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摘要 : ”电磁 环境 中 双 黑 洞 并 合 可 辐射 与 引力 波 同 频率 的 电磁 信号 ， 但 低频 电磁 波 有 可 能 被 星系 
介质 吸收 。 电 磁 波 在 穿越 大 气 过 程 中 会 被 吸收 ， 某 些 频 率 的 电磁 波 被 吸收 得 多 ， 某 些 频 率 的 电磁 
波 被 吸收 得 少 ， 被 吸收 得 少 的 那些 频率 窗口 被 称 为 电磁 波 的 大 气 窗 口 。 类 似 地 ， 电 磁 波 在 宇宙 空 
间 中 传播 会 被 星系 介质 吸收 。 着 重 探讨 电磁 波 在 银河 系 空间 中 传播 的 吸收 问题 。 研 究 发 现 ， 电 磁 
波 穿越 银河 系 空间 到 达 地 球 存在 下 限 截 止 频 率 。 天 球 上 不 同方 向 上 对 应 的 下 限 截止 频率 不 同 ， 存 
在 一 个 分 布 ， 该 分 布 与 银河 系 自由 电子 气 的 分 布 密切 相关 。 该 文通 过 对 该 截止 频率 分 布 进行 测定 
从 而 确定 银河 系 自由 电子 气 的 分 布 。 
关 键 词 : 电磁 波 ， 自由 电子 气 ， 截 止 频率 
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1 引 言 


引力 波 自 2015 年 9 月 被 美国 的 LIGO 引力 波 探 测 器 首次 探测 "到 以 来 ， 已 经 有 一 百 多 
例 引 力 波 事件 被 探测 到 ”。 以 引力 波 为 天 文 观测 手段 的 引力 波 天 文学 ”也 已 经 建立 起 来 ， 并 
以 惊人 的 速度 飞速 发 展 ”。 包 含 引 力 波 在 内 的 多 信使 天 文学 ”也 广 受 人 们 重视 。 以 2017 
年 8 月 观测 到 的 双 中 子 星 并 合 事 件 为 例 ， 引 力 波 与 电磁 波 的 联合 观测 同时 解决 了 短 伽 玛 身 
线 暴 的 中 心 引擎 问题 、 千 新 星 中 心 引 擎 问题 和 重 元 素 起 源 问题 三 大 天 文学 重大 科学 问题 。 这 
展示 了 引力 波 与 电磁 波 联合 观测 极 强 的 科学 发 现 能 力 。 

到 目前 为 止 ， 我 们 只 有 双 中 子 星 并 合 对 应 的 引力 波 与 电磁 波 联合 观测 结果 。 宇 宙 中 的 号 
洞 是 否 带 电 以 及 带 多 少 电量 是 近年 来 人 们 热烈 讨论 的 问题 一 ， 带 电 双 黑洞 并 合 “” 以 及 
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在 磁场 环境 下 的 双 黑 洞 并 合 “” 会 产生 与 引力 波 对 应 的 电磁 波 信号 ”。 值 得 注意 的 是 ， 这 
些 电磁 波 信号 具有 与 对 应 引力 波 信 号 相当 的 频率 。 以 LIGO 为 代表 的 地 面 引 力 波 探 测 器 对 
应 的 引力 波 频率 为 几 百 赫兹 ， 未 来 的 空间 引力 波 探测 器 对 应 的 引力 波 频率 为 之 赫兹 ， 脉 冲 星 
计时 对 应 的 引力 波 频 率 更 是 低 到 纳 赫兹 。 对 比 电磁 波 观测 ， 目 前 天 文 观测 使 用 的 电磁 波 最 低 
频率 也 在 10 MHz 以 上 。 典 型 的 射电 观测 在 GHz 频段。 

现 有 的 电磁 天 文 观测 没有 关注 10 MHz 以 下 频率 的 电磁 波 ， 其 主要 有 两 方面 的 原因 : 
(1) 10 MHz 以 下 频率 的 电磁 波 是 否 具 有 重要 和 有 趣 的 天 体 源 目前 尚 不 清楚 ; (2) 这 样 的 电磁 
波 能 否 穿 过 宇宙 空间 不 被 吸收 歼 尽 从 而 被 我 们 观测 到 也 是 关键 问题 。 根 据 我 们 前 面 的 描述 ， 
第 一 个 方面 不 再 是 问题 ， 目 前 我 们 具有 强烈 的 科学 动机 去 观测 10 MHz 以 下 其 低频 率 的 电磁 
波 信号 ， 第 二 个 方面 是 本 文 重点 关注 的 问题 。 

电磁 波 在 穿 过 地 球 大 气 的 过 程 会 被 吸收 是 我 们 已 经 非常 熟悉 的 物理 现象 ， 即 大 气 中 不 
同 分 子 对 不 同 频率 的 电磁 波 进行 一 定 程度 的 吸收 。 综 合 起 来 ， 整 个 大 气 对 不 同 频率 的 电磁 波 
产生 不 同 程度 的 吸收 效果 ， 因 此 形成 了 我 们 熟知 的 电磁 波 在 频率 空间 相对 大 气 的 透明 窗口 。 
即使 不 完全 透明 ， 我 们 也 对 不 同 频 率 的 电磁 波 的 吸收 系数 有 非常 全 面 的 了 解 。 

宇宙 空间 几乎 真空 ， 但 并 非 完 全 真空 ， 有 少量 的 星际 介质 存在 。 星 际 介质 的 主要 成 分 是 
自由 电子 气 ， 或 者 更 准确 地 说 ， 自 由 电子 气 和 离子 组 成 的 等 离子 体 。 等 离子 体 对 电磁 波 的 吸 
收 机 制 与 分 子 对 电磁 波 的 吸收 机 制 不 同 。 我 们 在 第 2 章 回顾 等 离子 体 同 电磁 波 的 相互 作用 
理论 和 吸收 机 制 。 在 第 3 章 我 们 把 该 吸收 机 制 应 用 到 银河 系 星际 介质 ， 分 析 计 算得 到 银河 
系 星际 介质 对 电磁 波 的 吸收 情况 ， 其 间 我 们 还 会 分 析 电 磁 波 吸收 情况 与 银河 系 自由 电子 气 
分 布 的 关系 。 最 后 在 第 4 章 ， 我 们 总 结 讨论 本 文 所 得 结果 ， 并 提出 通过 对 电磁 波 吸收 情况 
的 测定 来 确定 银河 系 自 由 电子 气 的 分 布设 想 。 


2 等 离子 体 与 电磁 波 的 相互 作用 


电磁 波 在 介质 中 传播 可 以 用 以 下 的 麦 殉 斯 书 方程 组 来 描述 : 


V (cE) = p ’ (1) 
V-B=0, (2) 
OB 
Vans ; (3) 
B O(c€E) 
Vx—= Ot ; (4) 


H 
其 中 ，p 是 自由 电荷 密度 ，7 是 传导 电流 密度 ，e 和 上 人 分 别 是 介质 的 介 电 系 数 和 磁 导 系数 ， 
而 瑟 和 B 分 别 是 电磁 波 的 电场 和 磁场 。 
对 于 自由 电子 气 ，e 和 jy 可 以 分 别 近 似 描述 四 为 
ulw) ~ ew) 、， wp 
a ~1, 7 = 1 -E (5) 
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a 


中 ，eo 和 pio 分 别 是 真空 介 电 系数 和 真空 磁 导 系数 ，w 是 电磁 波 的 频率 ，w EA 
气 对 应 的 等 离子 体 频率 : 


wp 一 人 /一 一 ， (6) 


其 中 , n 是 电子 数 密度 ，e 是 电子 电量 ，m 是 电子 质量 。 由 于 电子 数 密度 n 存在 空间 分 布 ， 
是 空间 位 置 的 函数 ， 所 以 对 于 星际 介质 u 基本 是 常数 ， 但 e 是 空间 位 置 的 函数 。 
考虑 到 常数 六 和 空间 函数 的 特点 ， 上 述 的 麦克 斯 韦 方 程 可 简化 为 : 


T -VE+V (£. e) =0. (7) 
对 于 单 频 电 磁 波 ， 
E = Etn” . (8) 
将 式 B) RAR 四， 我 们 得 到 变形 的 雍 姆 霍 兹 方程 : 
wpb +E +0 (EZ) <0 (9) 


星际 空间 自由 电子 气 密度 的 变化 很 慢 ， 变 化 尺度 在 0.001 pc 以 上 。 也 就 是 说 在 电磁 波 

源 和 探测 器 之 间 的 区 域 ， 如 果 我 们 把 区 域 分 成 很 多 尺度 小 于 这 个 变化 尺度 的 小 块 ， 在 每 个 小 
块 区 域内 ，n 可 以 近似 为 常数 。 于 是 式 ( 团 可 近似 为 通常 的 亥 姆 霍 兹 方程 : 

wepE +V°E=0 . (10) 


当 我们 把 上 述 方程 限制 到 我 们 关心 的 电磁 波源 所 在 方向 ， 并 用 x 表示 电磁 波 的 传播 方向 后 ， 
得 到 : 


weu + — =0 . (11) 

分 别 用 E; 和 Eii 代表 进入 第 i 块 小 区 域 前 的 电场 和 出 第 i 块 小 区 域 后 的 
第 i 十 1 块 小 区 域 前 的 电场 )。 我 们 可 以 通过 求解 式 (0) 得 到 关系 : 

En = Be ， (12) 

k=w/eu , (13) 


其 中 ，6z 是 第 i 块 小 区 域 的 x 厚度。 把 所 有 的 小 区 域 连 起 来 ， 我 们 得 到 探测 器 处 电场 Ba 
和 波源 处 电场 五 。 的 关系 : 


SN 
S 
~ 
Gt 
> 


Ey = E, aie kda ; (14) 


上 述 积分 从 波源 处 积分 到 探测 器 处 。 如 果 eS *a* 为 纯 实 数 ， 则 表明 中 间 没 有 吸收 ， 反 之 如 
Bot 的 虚 部 为 7， 则 表明 电磁 波 强度 由 于 吸收 衰减 e” 倍 。 
把 上 述 讨论 应 用 到 地 球 和 银河 系 边 缘 ， 我 们 就 得 到 电磁 波 被 银河 系 星 际 介质 的 吸收 
情况 。 
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3 银河 系 星 际 介质 对 电磁 波 的 吸收 


仿照 大 气 对 电磁 波 的 吸收 ， 我 们 定义 银河 系 星际 介质 对 电磁 波 的 吸收 系数 为 
E(w, 0,6) =1—e Un, (15) 


上 述 积分 沿 着 天 球 方向 (0,6) 从 地 球 积分 到 银河 系 边 缘 。 这 里 的 符号 S 代表 取 虚 部 。 本 文 
使 用 的 天 球 坐标 系 指 的 是 原点 在 太阳 ， 银 道 面 为 平面 xy， 太阳 银 心 连 线 为 轴 z 的 坐标 系 。 
这 个 坐标 系 不 是 惯性 坐标 系 ， 在 旋转 ， 但 旋转 速度 非常 慢 ， 大 约 2.2 亿 年 转 一 圈 。 

我 们 发 现 ， 吸 收 系数 El, 0, p) 相对 于 电磁 波 频 率 基 本 是 一 个 阶 越 函数 。 作 为 一 个 例子 ， 
我 们 任 取 一 个 方向 得 到 的 结果 (如 图 四 所 示 )， 也 就 是 说 在 给 定 的 方向 上 ， 存 在 一 个 电磁 波 
截止 频率 九 。 对 于 频率 大 于 该 截止 频率 的 电磁 波 不 会 被 银河 系 星际 介质 吸收 ， 反 之， 频率 小 
于 该 截止 频率 的 电磁 波 则 会 被 银河 系 星 际 介质 吸收 殉 尽 
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注 : 该 图 的 计算 基于 YMW16 电子 气 分 布 模型 。 


1 ”银河 系 星 际 介 质 对 电磁 波 的 吸收 系数 与 电磁 波 频率 的 关系 


银河 系 星际 介质 主要 是 自由 电子 气 。 很 多 学 者 对 银河 系 自由 电子 气 的 分 布 进行 了 仔细 
深入 的 研究 ， 代 表 性 结果 包括 Taylor 和 Cordes 提出 的 TC93 模型 号 。 在 过 去 对 于 自由 电 
子 分 布 轴 对 称 的 假设 基础 上 ，TC93 模型 研究 了 旋 臂 位 置 与 脉冲 星 分 布 之 间 的 相关 性 以 及 胶 
状 星云 色散 量 的 贡献 。TC93 模型 所 提供 的 已 知 脉冲 星 的 距离 估计 精度 约 为 25%。TC93 模 
型 在 银河 高 纬度 区 域 提 供 的 密度 偏 低 ， 只 可 推导 出 脉冲 星 距离 的 下 限 ， 且 在 旋 臂 内 部 的 电 
子 分 布 有 显著 的 误差 。 在 某 些 方向 上 提供 的 散射 太 少 ， 无 法 解释 脉冲 展 宽 。 建 立 在 TC93 BE 
型 基础 之 上 ，Cordes 和 Lazio 后 来 增加 了 对 更 多 测量 值 的 研究 并 且 改 进 了 对 其 星系 中 心 散 
射 强度 以 及 分 布 的 约束 ， 研 究 了 在 河 心 径 向 上 的 自由 电子 分 布 ， 研 究 了 局 部 悬臂 的 影响 。 将 
整个 银河 系 分 为 了 由 厚 盘 (thick disk), $; (thin disk) 和 旋 臂 (spiral arms) 三 个 主要 区 域 
分 别 进行 研究 。 在 计算 方法 上 使 用 了 迭代 似 然 法 (iterative likelihood method) 寻找 最 佳 拟 
合 模型 参数 。 由 此 ，Cordes 和 Lazio 7443] NE2001 模型 。 当 纬度 较 高 时 ，NE2001 模型 
的 准 ae 比 观 测 结果 偏 小 。 在 NE2001 模型 的 基础 上 ，Yao 等 人 ”将 整个 银河 系 分 
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成 厚 盘 、 薄 盘 、 旋 臂 、 星 系 中 心 (galactic center)、 局 域 泡 (local bubble region) 和 环 状 区 域 
(loop) 七 个 主要 区 域 ， 在 各 区 域 边 缘 时 研究 多 个 区 域 对 电子 密度 的 共同 影响 ， 提 出 YMW16 
模型 。YMW16 模型 在 95% 的 置信 区 间 范 围 内 优 于 NE2001 模型 近 40%. 

基于 YMW16 模型 ， 在 图 Ba) 我 们 给 出 电磁 波 截止 频率 在 不 同 天 球 方向 的 分 布 情况 。 
我 们 发 现 ， 银 河 系 星际 介质 吸收 电磁 波 对 应 的 截止 频率 在 5 x 103 ~ 1.1 x 104 Hz 之 间 。 
中 截止 频率 偏 高 的 方向 对 应 环 状 区 域 出 现 的 方向 。 作 为 对 比 ， 我 们 基于 NE2001 模型 ， 
Bib) 展示 了 电磁 波 截止 频率 在 不 同 天 球 方向 的 分 布 情况 。 对 比 图 欧 ) 和 图 b), RTE 
同 电子 气 模型 具有 明显 不 同 的 电磁 波 截止 频率 的 方向 分 布 情况 。 


截止 频率 /kHz 


a) 


11.86 203.68 


截止 频率 /kHz 
b) 


VE: a) 基于 YMW16 自由 电子 气 分 布 模型 ，b) 基于 NE2001 自由 电子 气 分 布 模型 。 


图 2 ”电磁 波 穿 过 银河 系 星 际 介质 到 达 地 球 的 截止 频率 在 天 球 上 的 分 布 
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系 而 到 达 地 球 ， 所 以 ， 我 们 不 可 能 观测 到 小 于 1000 Hz 频率 的 天 文 电磁 信号 
带电 黑洞 并 合 以 及 在 磁场 环境 中 并 合 双 黑 洞 的 引力 波 电磁 对 应 体 不 可 
成 的 射电 和 其 他 高 频 电 磁 信 号 男 当 别 论 。 
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由 于 银河 系 星际 物质 对 电磁 波 的 吸收 作用 ， 小 于 1000 Hz 频率 的 


月 


双 黑 洞 并 合 过 程 由 于 辐射 转移 等 
带电 或 磁场 环境 中 的 恒星 级 双 黑 洞 并 合 产 生 的 是 300 Hz ~ 3 kHz 的 特 低 


频 电磁 波 。 不 


由 于 银河 系 星际 介质 的 影响 ， 这 样 的 电磁 波 信和 号 无 法 传播 。 根 据 本 文 的 


磁场 环境 下 就 可 以 产生 3 x 103 ~ 3 x 104 Hz WER, Abe ee 


球 。 


=] 
“Jo 


而 超 导 量 子 干涉 器 件 (SQUID) 以 及 谐振 冷 线圈 等 技术 可 以 测量 到 这 样 的 


超 导 量 子 干涉 器 件 基 于 超 导 约瑟夫 森 效 用 和 磁 通 量子 化 效应 ， 对 这 个 频段 


电磁 波 信号 可 以 在 星际 空间 传播 。 如 果 存 在 质量 小 于 1 Mo 的 原初 黑洞 ， 


究 


际 介质 到 达 地 
微弱 其 低频 信 
的 磁场 分 量 非 


HRB. REA EAH SQUID 对 磁 分 量 的 测量 已 达到 10- 工 水 平 。 由 于 其 极 高 的 灵敏 


技术 研究 所 已 经 开始 精密 


并 绘制 电磁 波 
认 银 河 系 自 由 


FE, SQUID 已 经 用 于 一 些 空间 物理 实验 中 ， 如 Gravity Probe B™ 和 暗物质 的 探测 。 
目前 ， 利 用 超 导 量 子 干涉 器 件 ， 中 国 科 学 院 上 海 微 系统 与 信息 

测量 地 球 大 气 层 环境 中 的 其 低频 信号 ” 。 我 们 期 待 未 来 能 利用 SQUID 测量 
截止 频率 在 不 同 天 球 方向 的 分 布 情况 。 这 样 的 测量 结果 可 以 很 好 地 区 分 和 证 
电子 气 分 布 模型 。 
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Very Low Frequency Electromagnetic Wavse Absorption and 
Free Electron Gas Distribution in the Milky Way 
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Abstract: Electromagnetic waves will be absorbed when they pass through the atmo- 
sphere. Electromagnetic waves with some frequencies are absorbed more, while electromag- 
netic waves with other frequencies are absorbed less. The frequency range within which 
the electromagnetic waves are less absorbed is called the atmospheric window of electro- 
magnetic waves. Similarly, the propagation of electromagnetic waves in the universe will be 
absorbed by the galactic medium. This paper investigates the problem of the absorption 
of electromagnetic wave when it propagates in the Milky Way. Our analysis shows that 
there is a lower cut-off frequency for the electromagnetic waves which can reach the earth 
passing through the Milky Way. Interestingly, the lower cut-off frequency is different respect 
to different direction on the celestial sphere, which gives a distribution. This distribution 
is closely related to the distribution of free electron gas in the Milky Way. We propose to 
determine the distribution of free electron gas in the Milky Way by measuring the cut-off 


frequency distribution. 


Key words: electromagnetic wave; free electron gas; lower cut-off frequency 


